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^— N ' Summary 

(N ■ 

f^ . A universal initial mass function (IMF) is not intuitive, but so far no convincing evidence for 

~y \ a variable IMF exists. The detection of systematic variations of the IMF with star- forming 
conditions would be the Rosetta Stone for star formation. 
^^ . In this contribution an average or Galactic-field IMF is defined, stressing that there is 

evidence for a change in the power-law index at only two masses: near 0.5 M0 and 0.08 Mq. 

04 ' Using this supposed universal IMF, the uncertainty inherent to any observational estimate of 

J^ , the IMF is investigated, by studying the scatter introduced by Poisson noise and the dynamical 

f— ^ ' evolution of star clusters. It is found that this apparent scatter reproduces quite well the 

^^ ' observed scatter in power-law index determinations, thus defining the fundamental limit within 

Q^ \ which any true variation becomes undetectable. The absence of evidence for a variable IMF 

^— V ' means that any true variation of the IMF in well studied populations must be smaller than this 

^^ . scatter. 

fH \ Determinations of the power-law indices a are subject to systematic errors arising mostly 

f^ 1' from unresolved binaries. The systematic bias is quantified here, with the result that the 

Q . single-star IMFs for young star-clusters are systematically steeper by Aa k, 0.5 between 

H \ 0.1 and IMq than the Galactic-field IMF, which is populated by, on average, about 5 Gyr 

^ ■ old stars. The MFs in globular clusters appear to be, on average, systematically flatter than 

. . ! the Galactic- field IMF (Piotto & Zoccah 1999; Paresce & De Marchi 2000), and the recent 

. 5^ I detection of ancient white-dwarf candidates in the Galactic halo and absence of associated 

j^ ' low-mass stars (Mendez & Minniti 2000; Ibata et al. 2000) suggests a radically different IMF 

3 , for this ancient population. Star-formation in higher-metallicity environments thus appears to 

■ " " ' produce relatively more low-mass stars. While still tentative, this is an interesting trend, being 
consistent with a systematic variation of the IMF as expected from theoretical arguments. 



Subject headings: stars: mass function - stars: formation - binaries: general - open clusters and 
associations: general - globular clusters: general - stellar dynamics 



1. INTRODUCTION 

Fundamental arguments suggest that the initial mass function (IMF) should vary with the pressure and 
temperature of the star-forming cloud in such a way that higher-temperature regions ought to produce 
higher average stellar masses (Larson 1998). This is particularly relevant to the formation of population III 
stars, because the absence of metals and more intense ambient radiation field imply higher temperatures. 



The IMF inferred from Galactic-field star-counts can be conveniently described by a 3-4 part power-law 
(eqs. |l| and H below). The Galactic field was populated by many different star-formation events. Given 
this well-mixed nature of the solar neighbourhood, present-day star-formation ought to lead to variations 
about the Galactic- field IMF. In particular, a systematic difference ought to be evident between low-density 
environments (e.g. Taurus- Auriga; p Oph) and high-density regions (e.g. Orion Nebula Cluster, ONC), 
because above a certain critical density, star-forming clumps interact with each other before their collapse 
completes (Murray & Lin 1996; Allen & Bastien 1995; Price & Podsiadlowski 1995; Klessen & Burkert 
2000). On considering the ratio between the fragment collapse time and the collision time-scale and applying 
the analysis of Bastien (1981, his eq. 8), it becomes apparent that the IMF in clusters similar to p Oph 
cannot be shaped predominantly through collisions between collapsing clumps. This is supported through 
the finding by Motte, Andre & Neri (1998) that the pre-stellar-clump MF in p Oph is similar to the observed 
MF for pre-main sequence stars in p Oph. It is somewhat steeper than the Galactic-field IMF, especially 
if the binary systems that must be forming in the pre-stellar cores are taken into account. Noteworthy is 
that both, the pre-stellar clump MF and the Galactic-field IMF have a reduction of the power-law index 
below about 0.5 Af0. In the core of the ONC, however, pre-stellar cores most likely did interact significantly 
(Bonnell, Bate & Zinnecker 1998; Klessen 2001). Furthermore, once the OB stars ignite in a cluster such as 
the ONC, they have a seriously destructive effect through the UV flux, strong winds and powerful outflows, 
and so are likely to affect the formation of the least massive objects, including planets. This can happen for 
example through destruction of the accretion envelope, so that extreme environments like the Trapezium 
Cluster may form a surplus of unfinished stars (brown dwarfs, BDs) over Taurus-Auriga. Luhman (2000) 
finds empirical evidence for this, but detailed dynamical modelling is required to exclude the possibility 
raised here that at least part of this difference may be due to the disruption of BD-BD and star-BD 
binaries in a dynamically evolved population such as the Trapezium Cluster. 

A conclusive difference has not been found between the IMF in Taurus- Auriga (Kenyon & Hartmann 
1995; Briceno et al. 1998) and p Oph (Luhman & Rieke 1999) on the one hand, and the ONC (Palla & 
Stabler 1999; Muench, Lada & Lada 2000; Hillenbrand & Carpenter 2000) on the other. Similarly, Luhman 
& Rieke (1998) point out that no significant IMF differences for pre-main sequence populations spanning 
two orders of magnitude in density have been found. Such conclusions rely on pre-main sequence tracks 
that are unreliable for ages younger than about 1 Myr (I. Baraffe, priv.commun.), because the density, 
temperature and angular momentum distribution within the pre-main sequence star is likely to remember 
the accretion history (Wuchterl & Tscharnuter 2000). Nevertheless, in support of the universal-IMF notion, 
it is remarkable how similar the Galactic field MF is to the MF inferred for the Galactic bulge (Holtzman et 
al. 1998; Zoccali et al. 2000), again with a fiattening around 0.5 Mq. Presumably star-formation conditions 
during bulge formation were quite different to the conditions witnessed in the Galactic disk, but the bulge 
and disk metallicities are similar. Further related discussions on this topic can be found in Gilmore & 
Howell (1998). 

The quest for detecting variations in the IMF has been significantly pushed forward by Scalo (1998), 
who compiled determinations of the logarithmic power-law index, F (eq. y), for many clusters and OB 
associations in the Milky Way (MW) and the Large Magellanic Cloud (LMC), which has about l/5th— l/3rd 
the metallicity of the MW (e.g. Holtzman et al. 1997). While no systematic variation is detectable in a plot 
of r against stellar mass, m, between populations belonging to the two galaxies, a large constant scatter in 
F for stars more massive than 1 Mq is evident instead. This raises the question of how large apparent IMF 
variations are due to small number statistics and other as yet unexplored observational uncertainties, and 
if this noise can mask, or even render undetectable, any true variations of the IMF. 



Elmegreen (1999) shows that statistical variations of T, that are not dissimilar to the observed ones, 
result naturally from a model in which the Salpeter IMF constructs from random sampling of hierarchically 
structured clouds, if about N = 100 stars are observed. This model predicts that the scatter in F must 
decrease with increasing N. 

In this contribution the reductionist philosophy is followed according to which all non-star-formation 
sources of apparent variations of the IMF must be understood before the spread of F can be interpreted 
as being due to the star formation process. To achieve this, an invariant IMF is assumed to study three 
possible contributions to the large scatter seen in the alpha-plot: (i) Poisson scatter due to the finite 
number of stars in a sample. This is similar to Elmegreen's approach, except that no explicit link to the 
distribution of gas clumps is made, (ii) Loss of stars of a preferred mass-scale as their parent star-clusters 
evolve dynamically. This dynamical loss is not a simple function of stellar mass, because of the complex 
stellar-dynamical events in a young cluster. For example, while low-mass members preferably diffuse 
outwards as a result of energy equipartition, massive stars sink inwards where they meet and expel each 
other rather effectively. Finally, (iiia) wrong mass estimates because most stars are born in binary systems, 
and observers usually cannot resolve the systems. The simplest approach, taken here, is to replace the two 
component masses by the combined mass of the binary system, and to measure the system MF. 

Issues also contributing to the scatter but not dealt with here are the following: (iiib) An observer 
infers the mass of a star from the observed luminosity incorrectly if the star is an unresolved binary, (iiic) 
wrong mass estimates from luminosities in the event of higher-order multiplicities, which is a major bias 
for massive stars (e.g. Preibisch et al. 1999), and (iva) stellar evolution and the application of incorrect 
pre-main sequence and main-sequence evolutionary tracks, which corrupts the masses inferred from observed 
quantities as the stars evolve to or along the main sequence, and (ivb) incorrect estimates of stellar masses 
as a result of rapidly rotating massive stars and the use of non-rotating stellar evolution models. One issue 
to be stressed in this connection is that stellar evolution theory retains significant uncertainties (Kuruzc 
2000; Maeder & Meynet 2000; Heger & Langer 2000), which can only be reduced through continued 
attention. 

The present study thus probably underestimates the scatter by focusing on points (i) to (iiia), but 
allows an assessment of the fundamental limits within which apparent IMF variations swamp true variations. 

The alpha-plot and the form of the universal IMF adopted here are introduced in Section 0, and 
statistical variations of the power-law index are studied in Section ^. The star-cluster models are described 
in Section ^ and Section ^ contains the results on the variation of the MF. In Section |6|the dichotomy in 
the alpha-plot and available evidence for a truly variable IMF are discussed. The conclusions are presented 
in Section j^. 



2. THE ALPHA-PLOT AND THE GALACTIC-FIELD IMF 

Observational data in the alpha-plot are used to infer a universal IMF. 

2.1. The alpha plot 

Scalo (1998) combined available IMF estimates for star clusters and associations by plotting the 
power-law index, F (eq. below), against the mean logiorn of the mass range over which the index is 



measured (his fig. 5). Fig. m shows these same data by plotting the power-law index a = 1 — F against 
logioTOav The alpha-plot clearly shows the flattening of the IMF for ?7i ^ 0.5 Mq. It also shows no systematic 
difference between MW and LMC populations, as already shown by Massey et al. (1995b) for massive 
stars. This is also verified for 0.6^rn < 1.1 A/q by Holtzman et al. (1997), who use deep HST photometry 
for LMC fields and apply Monte Carlo models that include binary systems, various star-formation histories 
(sfh) and metallicities as well as observational errors. 

The models discussed in Section g show that unresolved binary systems mostly affect the region 
TO ^ 1 Mq , the data for which are listed in Table |^. Perusal of the references shows that only Meusinger et 
al. (1996) attempted a correction for unresolved binary systems. However, they adopted an artificial model 
of Reid (1991; see discussion in Kroupa 1995a), in which the binary proportion is only 40 per cent, half of 
which have similar companion masses. This is an unrealistic model (Kahler 1999), and leads to essentially 
no difference between the system and single-star LFs (fig. 7 in Meusinger et al. 1996, compare with fig. 11 
in Kroupa 1995d). Their binary correction can thus be safely neglected. 



2.1.1. m>3MQ 

For m<^3A/0 the data suggest that the Salpeter power-law value, a = 2.3, is a reasonable fit over 
the whole range, as is also stressed by Massey (1998). Massey & Hunter (1998), for example, deduce 
that a « 2.3 for 2.8 < m < 120 Mq in the massive cluster R136 in the LMC. This value is thus adopted 
throughout the rest of this paper, although notable examples of exotic clusters exist. The two massive, 
apparently young (2-4 Myr) Arches and Quintuplet clusters lying very close to the Galactic centre 
(projected distance 30 pc) have a « 1.65 (Figer et al. 1999), and the Galactic star-burst template cluster 
NGC 3603 is found to have a « 1.7 (Eisenhauer et al. 1998). Further work is desired to establish the exact 
nature of the central clusters, and clarify the age discrepancy between the low-mass and massive stars noted 
for NGC 3603, a problem possibly associated with pre-main sequence models. 

It is important to keep in nnnd that a may be systematically steeper than a = 2.3 (or 1.7) due 
to unresolved binary systems, which are not usually corrected for in IMF estimates. The multiplicity 
proportion of massive stars is very high (Mason et al. 1998). For example, Preibisch et al. (1999) find 
that the OB stars in the well studied ONC have, on average, 1.5 companions. For each primary, there 
is thus usually more than one secondary that adds at lower masses, steepening the observed IMF when 
corrected for. The effect depends on a, and Sagar & Richtler (1991) calculate that Aa = -1-0.34 for a = 2.5 
and a binary proportion / = 0.5 (eq. ra below). If / = 1 (each massive primary has 1.0 companions) they 



logiomav a cluster ref. 

(Williams et al. 1995b) 
(Comeron, Rieke & Rieke 1996) 
(Comeron, Rieke & Rieke 1996) 
(Williams, Rieke & Staufler 1995a) 
(Meusinger et al. 1996) 
(Hillenbrand 1997) 
(Williams, Rieke & Staufler 1995a) 

Table 1: The data from Scale's (1998) compilation with rriav < 1 Mq 
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Fig. 1. — The alpha-plot. The data show the compilation by Scalo (1998) of determinations of a over different 
mass ranges for Milky- Way (MW) and Large-Magellanic-Cloud (LMC) clusters and OB associations. Unresolved 
multiple systems are not corrected for. The large open triangles (Muench, Lada & Lada 2000 from Orion Nebula 
Cluster observations, binary corrections not applied) serve to illustrate the present knowledge for m < 0.1 M©. 
The horizontal long-dashed lines in the BD regime are the Galactic-field IMF (eq. Gl) with associated approximate 
uncertainties. For 0.08 < m < 1.0 M© the thick short-dashed lines represent the KTG93 single-star IMF (Kroupa, 
Tout & Gilmore 1993), which has 0:3 = 2.7 for m > 1 Mq from Scalo's (1986) determination. The long-dashed lines 
for m > I Mq show the approximate average a — 2.3, which is adopted in the Galactic-field IMF (eq. 0). The Miller 
& Scalo (1979) log-normal IMF for a constant star-formation rate and a Galactic disk age of 12 Gyr is plotted as 
the diagonal long-dash-dotted line. The long-dash-dotted horizontal lines labelled "SN" are those 03 — 0.70(1.4) for 
which 50 % of the stellar (including BD) mass is in stars with 8 — 50(8 — 120) Mq. The vertical dotted lines delineate 
the four mass ranges (eq. 0), and the shaded areas highlight those stellar mass regions where the derivation of the 
IMF is additionally complicated due to unknown ages, especially for Galactic field stars: for 0.08 < m < 0.15 M© 
long-contraction times make the conversion from an empirical LF to an IMF dependent on the precise knowledge of 
the age, and for 0.8 < m < 2.5 Mq post-main sequence evolution makes derived masses uncertain in the absence of 
precise age knowledge. A few of the MW data are labelled by their star-clusters, and Table 111 lists the niav < 1 Mq 
data. 



obtain Aa = +0.40. Aa is likely to be larger still, because each massive primary probably has more than 
one companion, typically. Since the empirical data in Fig. |l| implies an average a ~ 2.3 for m^3 Mq, the 
true single-star IMF may in fact have a « 2.7(= 2.3 + 0.4), or even larger. A similar conclusion is reached 
by Scalo (1998, at the end of his section 4). Such corrections will not be removed if spectroscopic mass 
determinations are used instead of the inferior mass-estimates using photometry (Massey et al. 1995a), 
since unresolved systems will have similar effects on a spectroscopic sample. 

In this paper the approximate average a = 2.3 is adopted, with the aim of studying the effect of 
unresolved binary systems on the a inferred from the system MF, which an observer would deduce from the 
mixture of single stars and binary systems in a population resulting from star-cluster evolution with initially 
/ = 1 . Because the assumptions (Section |4.3| ) imply that massive stars have very-low mass companions in 
this model, and because only binary systems are searched for in the data reduction software, the resulting 
model bias will be an underestimate. Further work is necessary to address this particular issue, which is 
also discussed further at the end of Section |^ and in Section iA . 



The remarkable feature for to ^ 3 Mq in the alpha plot is the constant scatter, and that the various 
power-law indices are distributed more or less randomly throughout the region a = 2.3 ± 0.7 without a 
significant concentration towards some value. 



2.1.2. 0.8<to<3Mq 

The region for 0.8 ^m^3 Mq shows an unusually large scatter. It is shaded because this particular 
mass range is problematical for a number of reasons. 

Analysis of Galactic-field star-counts run into the difficulty that the age of the Galactic disk is 
comparable to the life-time of these stars so that stellar evolution corrections become very significant, but 
for this the sfh must be known (Scalo 1986; Haywood, Robin & Creze 1997). That interesting constraints 
can be placed on the MW IMF using independently derived sfhs is shown by Maciel & Rocha-Pinto (1998), 
where the problems associated with the estimation of the field- IMF for massive stars are documented. 

The large spread of the cluster values in the region 0.8 ^ to, ^ 3 Mq may be due to the fact that the 
observed clusters have ages such that the stars in this mass range count to the most massive remaining 
in the clusters. They are thus subject to advanced stellar evolution and/or dynamical ejection from the 
cluster, because the most massive stars usually interact in the vicinity of the cluster core. Which of these 
is applicable is a sensitive function of the age of the cluster and the number of stars in it (more on this 



in Section 5_^). Finally, stellar evolution is by no means a solved subject for stars in this mass range 
(Dominguez et al. 1999) with remaining significant uncertainties. This compromises the conversion of 
stellar luminosity to mass. Ignoring the large scatter in this mass range, it can be seen that a single 
power-law index becomes applicable for ttt, > 0.5 Mq. 



2.1.3. 0.1<TO<1M, 







The Galactic- field single-star IMF fits the data shown in Fig. ^ exceedingly well for 0.1 < to < IMq 
(that this agreement may be fortuitous though is shown in Section |6.2D . In particular, it is remarkable that 
the data suggest a change in a near 0.5 Mq, as was initially derived from solar-neighbourhood star-counts 
by Kroupa, Tout & Gilmore (1991, hereinafter KTG91), and later confirmed by Kroupa, Tout & Gilmore 



(1993, hereinafter KTG93) and Kroupa (1995a) using a different mass-luminosity relation, a much more 
detailed star-count analysis including main-sequence and pre-main sequence stellar evolution, and with 
different statistical tests. Similar work by Gould, Bahcall & Flynn (1997) using HST star-counts and Reid 
& Gizis (1997), who study a proposed extension of the nearby stellar sample to somewhat larger distances, 
also confirm these findings, as do Chabrier & Baraffe (2000), who estimate a « 1.2 ± 0.1 using the nearby 
volume-limited LF. 

Of special importance is the mass range 0.5 — 1 Mq: The local sample of known stars is sufficiently 
large in this mass range that the nearby volume-limited LF is very well defined (Kroupa 2001a). Also, 
unresolved binaries do not significantly affect the LF in this mass range, because the stellar sample does 
not contain stars with to > 1 Mq that can hide lower-mass companions. The mass-luminosity relation is 
also well understood for these stars, so that the MF determination should be accurate and precise. It is 
not surprising that the power-law slope has changed little over the decades (Salpeter 1955: a = 2.35 for 
0.4 < m < 10 Mq). From Fig. |l] an uncertainty of a = 2.3 ± 0.3 is adopted. 

Unfortunately, the local sample of stars with to,^O.5M0 is incomplete for distances larger than 
d w 5 pc, in contradiction to the belief by Reid & Gizis (1997), who use spectroscopic parallax measurements 
to extend their proposed volume-limited sample using previously- known stars. Malmquist bias pollutes their 
sample by multiple systems that are much further away. The seriousness of the incompleteness of the nearby 
stellar census is shown by Henry et al. (1997), and is also pointed out by Chabrier & Baraffe (2000). This 
situation can only be improved with large-scale and deep surveys that find candidate nearby M dwarfs with 
subsequent trigonometric parallax measurements to affirm the distance, such that a volume limited sample 
can be constructed. This will be possible through the upcoming astrometric space missions DIVA (Roser 
1999) and GAIA (Gilmore et al. 1998). Being aware of this situation, the KTG studies combined the local 
{d < 5.2 pc) volume limited sample with flux limited deep photometric surveys, performing detailed Monte 
Carlo modelling of both Galactic-field samples. This pedantic separation of the two star-count samples 
is necessary as completely different biases and errors operate. The result is the conservative uncertainty 
range of a = 1.3 ± 0.5 for 0.08 - 0.5 Mq (KTG93). That the Galactic-bulge MF shows an indistinguishable 
behaviour to the Galactic-field MF in this mass range was already pointed out in Section |]. 



2.1.4. TO < 0.08 Mq 

For sub-stellar masses the constraints have improved dramatically in the past few years as a result of 
the significant observational effort and instrumental advances. In the ONC, Muench, Lada & Lada (2000) 
and Hillenbrand & Carpenter (2000) find — 1 ^a^ 1, although the pre-main sequence tracks are unreliable 
at these ages. Similarly, in p Oph Luhman & Rieke (1999) estimate a « 0.5, which is also found for IC348 
by Najita, Tiede & Carr (2000). In the Pleiades Cluster, Martin et al. (2000) estimate a « 0.53, and for 
the solar neighbourhood, Reid et al. (1999) quote 1 ^a^2, whereas Herbst et al. (1999) estimate a < 0.8 
with 90 per cent confidence on the basis of no detections but accounting correctly for Galactic structure. 
For the time-being, a — 0.3 ± 0.7 is a reasonable description of the IMF for BDs, and it will be shown in 



Section 5.2 that the observed MF depends sensitively on the dynamical age of the population. 

The region 0.08 — 0.15 Mq is shaded in Fig. |l| to emphasise the uncertainties plaguing Galactic-field 
star-count interpretations as a result of the long pre-main sequence contraction times for these stars. 
As with the 0.8 — 3Mq region, the sfh must be known. The sfh has most recently been constrained by 
Rocha-Pinto et al (2000). 



2.2. The universal IMF 

The available constraints can be conveniently summarised by the multiple-part power-law IMF (see 
Kroupa 2001b for details), 

^(m) ex TO~"* = m'^^ , (1) 

where 

ao = +0.3 ± 0.7 , 0.01 < m/MQ < 0.08, 

ai=+1.3±0.5 , 0.08 < m/Me < 0.50, 

as = +2.3 ±0.3 , 0.50 < m/M© < 1.00, ^' 

as = +2.3 ±0.7 , l.OO<m/M0, 

and £,{m) dm is the number of single stars in the mass interval m to to + dm. The uncertainties correspond 
approximately to 99 per cent confidence intervals for m^O.5M0 (Fig. nl), and to a 95 per cent confidence 
interval for O.f — 0.5 M© (KTG93). Below 0.08 M© the confidence range is not well determined. 

Note that this form differs from Scalo's (1998) recommendation, mostly because the correct structure 
in the luminosity function below I Mq is accounted for here. There is evidence for only two changes in 
the power-law index, namely near 0.5 Mq and near 0.08 Mq. The frequently used Miller & Scalo (1979) 
IMF fails in the region 0.5 — IMq, and especially for m^bMQ (Fig. nl see also Fig. RJ below). A useful 
representation of the IMF is achieved via the logarithmic form, 

Cl("i) =^(m)lnlOTO, (3) 

where ^l dlog^Qm ex mF' dlog^Qm = m~^' dlogiQm is the number of stars in the logarithmic mass interval 
logioTO to logiom + dlogjQTO. 

The adopted IMF (eq. g) has a mean stellar mass < m >= 0.36 Mq for stars with 0.01 < m < 50 Mq, 
and leads to the following stellar population: 37 % BDs (0.01 — 0.08 M0) contributing 4.3 % to the stellar 
mass, 48 % M dwarfs (0.08-0.5 Afg) contributing 28 % mass, 8.9 % "K" dwarfs (0.5- 1.0 M©) contributing 
17 % mass, 5.7 % "intermediate mass (IM) stars" (1.0 - 8.0 Af©) contributing 34 % mass, and 0.37 % "O" 
stars (> 8 Mq) contributing 17 % mass. 

A remarkable property of eq. g is that 50 per cent of the mass is in stars with 0.01 < to < 1 Mq and 
50 per cent in stars with 1 — 50 Af©. Also, if a4 = 0.70 (m > 8 Af©) then 50 per cent of the mass is in stars 
with 8 < TO < 5OAf0, whereas a4 = 1.4 implies 50 per cent mass in 8 — 120 A/© stars. These numbers are 
useful for the evolution of star clusters, because supernovae (SN) lead to rapid mass loss which can unbind 
a cluster if too much mass resides in the SN precursors. This is the case in clusters that have a^ = 1.80: 
stars with 8 < tti < 120 Af© contain 53 per cent of the mass in the stellar population! It is interesting that 
a w 1.8 for TO^ 1 Af© forms the lower bound on the empirical data in Fig. nl But even 'normal' {a — 2.3) 
star clusters suffer seriously through the evolution of their to > 1 Af© stars (de La Fuente Marcos 1997). 



3. PROCEDURE AND STATISTICAL VARIATION 

One contribution to the scatter seen in the alpha-plot (Fig. H) is Poisson noise. This can be studied by 
sampling N stars from the adopted IMF (eq. ||), and studying the variation of a with N. 

In order to construct synthetic alpha-plots, the following procedure is adopted. N masses are obtained 
by randomly sampling eq. with lower mass limit m\ = 0.01 Af© and upper mass limit m^ — 50 Af©. 
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This upper mass limit is chosen for consistency with the stehar-dynamical models (Section ^). The MF is 
constructed by binning the masses, m, into 30 logio?7i bins which subdivide the range —2.1 < logj^o'^ — +2.1. 
Power-laws are fitted using weighted linear regression (e.g. Press et al. 1994) to sub-ranges that are defined 
as follows: 

logio(O.Ol) < Im < logio(0.08), 

logio(0.08) <\m< logio(0.50), 

logio(0.40) < l™ < logio(1.20), 

logio(l.OO) < l™ < logio(3.50), ^ ' 

logio(3.00) < Im < logio(9.00), 

logio(5.00) < Im, 

where Im = \ogiQm/MQ, and the numbers, nb, behind the mass-range number (e.g. 61) are the number of 
mass bins in the histogram in that particular mass range (e.g. nbi — 6). This sub-division ensures that the 
different mass regions in which a^ is known to be constant (eq. 0) are not mixed up, but also allows studying 
the fitted a at values of Im where, for example, stellar evolution and/or dynamical effects are expected to 
be important. The result is a(l?7iav), where l?7iav is the average of Im over the nbj (j = 1, 6) bins. In cases 
where the number of stars is too small, or the highest mass star is less massive than 10^'^ Mq, some of the 
highest-mass bins remain empty, causing l?Tiav in mass-range 66 to vary between renditions. 

The IMF is plotted together with two renditions using N = 10"^ stars in Fig. g, to illustrate the 
procedure. The resulting alpha-plot is shown in Fig. y| for many more renditions and different N. The 
input IMF is obtained essentially exactly for N = 10^ and 10^, verifying the procedure. The figure shows 
that deviations begin to occur for N = 10'* in the two highest mass-ranges (65 and 66), because these 
contain only a few per cent of N , i.e. a few hundred stars, spread over about 10 mass bins. For smaller 
N the scatter of a(lTOav) becomes larger, with the average reproducing the IMF except when the MF is 
under-sampled at large masses. 

Fig. H illustrates this sampling bias. The under-sampling of the histogram in the highest mass bins, 
when iV is too small, leads to an apparent flattening of the MF in the most massive bins accessible to the 
stellar population, as is evident in Fig. |[ It is also evident in fig. 2 of Elmegrecn (1999), and in typical 
star-count data, such as used by Massey et al. (1995b, their fig. 5) to infer the power-law index. Such 
samples contain typically a few dozen stars only (their table 5). This is interesting, possibly implying that 
the correct single-star IMF may be steepening, i.e. have an increasing a, with \m at the largest masses, 
since the uncorrected data suggest a constant a for m ^ 1 Mq. This issue, together with the bias through 
the high multiplicity fraction, will require more explicit modelling of the biases affecting the observed IMF 
for massive stars. 

In conclusion, Fig. ^ demonstrates that the observed scatter is arrived at approximately for populations 
that contain 10^ ^N ^ 10^ stars, which is quite typical for the type of samples available. 



4. STAR CLUSTER MODELS 

In Section g apparent variations of the IMF are discussed that result purely from statistical fluctuations. 
Additional sources of uncertainty are listed in the Introduction. Section |5| concentrates on quantifying the 
apparent variations that arise from stellar-dynamical effects and unresolved binary systems. To achieve 
this, a range of star-cluster models are constructed. This approach is relevant to populations in young 
clusters, OB associations and even the Galactic field, because most stars form in embedded clusters (Lada 
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logio(m) [Mq] 

Fig. 2. — The adopted logarithmic IMF (eqs. g and |3|), ^l/10^, for 10^ stars (solid histogram). Two random 
renditions of this IMF with 10 stars are shown as the heavy and thin dotted histograms. The mass-ranges over 
which power- law functions are fitted are indicated by the arrowed six regions (eq. m, while thin vertical dotted lines 
indicate the masses at which at changes. 




logi(,m [Mq] 



Fig. 3. — Purely statistical variation of a in the six mass ranges (eq. H) for different A'^ as indicated in the key. 
Large outer squares indicate those a fits obtained with nb — 2 and 3 mass bins. The open circles, open triangles, 
vertical and horizontal lines are as in Fig. nl. 
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& Lada 1991; Kroupa 1995b). 



4.1. Codes 

The dynamical evolution of the clusters studied here is calculated using Nbody6 (Aarseth 1999; 
Aarseth 2000), which includes state-of-the art stellar evolution (Hurley, Pols & Tout 2000), a standard 
Galactic tidal field (Terlevich 1987), and additional routines for initiating the binary-rich population 
(Kroupa 1995c). The A^-body data are analysed with a large data-reduction programme that calculates, 
among many quantities, the binary proportion and MFs. 



4.2. The clusters 

The cluster models are set-up to have the same central density, pc = 10^'* stars/pc"^, as observed 
in the Trapezium Cluster (McCaughrean & Stauffer 1994), giving a half-mass diameter crossing time 
Across — 0.24 Myr. The centre of masses of the binary systems follow a Plummer density distribution 
(Aarseth, Henon & Wielen 1974) with half- mass radius i?o.5- The average stellar mass is independent of the 
radial distance, i?, from the cluster centre, and the clusters are in virial equilibrium. Their parameters are 
listed in Table |[ Cluster evolution is followed for 150 Myr. 



4.3. The stellar population 

Stellar masses are distributed according to the IMF (eq. g) with m\ — 0.01 A/q and m^ — 50 Mq. This 
upper mass limit is half as large as the upper limit on the mass range used to evaluate the MP (eq. |j) , to 
take into account stellar mergers. Merging can occur during pre-main sequence eigenevolution, as detailed 
below. The default models assume as = a2, but one model is also constructed with the possibly more 
realistic value a^ = 2.7 (this model has ao = 0.5 for historical reasons). 



model 


iV 


iVbin 


Ro.5 


<m> 


0"3D 


ircl 


K, 








[pc] 


[Mq] 


[km/s] 


[Myr] 




B800 


800 


400 


0.19 


0.4 


1.6 


0.8-1.4 


5 


B3000 


3000 


1500 


0.30 


0.4 


2.5 


2.4-4.4 


5 


B1E4 


10* 


5000 


0.45 


0.4 


3.7 


6.8-12.5 


2 


BlE4d 


10* 


5000 


0.45 


0.3 


3.2 


7.9-14.5 


2 



Table 2: Cluster models: A'^ and A'^bin are the initial number of stars and binaries in each model (not taking into 
account mergers), Ro.s is the half-mass radius, and <m> is the average stellar mass. The three-dimensional velocity 
dispersion is (T3d, the median relaxation-time is irci. It's range results from assuming / = 1 and / = 0, respectively, 
since / evolves. The number of calculations per model is A'run. Model BlE4d has as = 2.7 (m > IMq, eq. bl), 
whereas the other (default) models have as — a2 — 2.3. It took about 4 months to assemble these data on standard 
desk-top computers. 
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Binaries are created by pairing the stars randomly. The binary proportion is 

where iVsing and iVbin are the number of single-star and binary systems, respectively. A birth binary 
proportion / = 1 is assumed. The initial mean system mass is 2 <m>, with <m> being the average 
stellar mass. This results in an approximately flat mass-ratio distribution (fig. 12 in Kroupa 1995c). Note 
however that encounters in clusters lead to the preferred disruption of binaries with low-mass companions. 
The initially "random" mass-ratio distribution evolves rapidly towards a distribution in which low-mass 
companions are less frequent, but still preferred (Kroupa 1995c). This is consistent with observations in 
that G-dwarf primaries (Duquennoy & Mayor 1991), Cepheids (4 — 9Mq, Evans 1995) and possibly OB 
stars (Mason et al. 1998; Preibisch et al. 1999) prefer low-mass companions. 

Periods and eccentricities are distributed following Kroupa (1995c). The periods range from about 1 d 
to 10^ d, and pre-main sequence eigenevolution changes the periods, mass ratios and eccentricities such 
that they are consistent with observational constraints for late-type main sequence stars with short periods. 
Eigenevolution is the collective name for system-internal processes that evolve the orbital parameters, such 
as tidal circularisation, mass transfer, and interactions with circum-stellar and circum-binary discs. One 
feature of the pre-main sequence eigenevolution model is that secondary companions gain mass during 
accretion if the peri-astron distance is smaller than a critical value. This affects the IMF by slightly 
reducing the number of low-mass stars, and slightly increasing the number of massive stars. Also, in some 
rare cases the binary companions merge giving iVging > 0, so that the true initial binary proportion is less 
than unity. Since only short-period binaries are affected by eigenevolution, the overall changes to the IMF 
are not significant. 

The resulting single-star and system MFs are shown in Fig. Q. This figure demonstrates that the IMF 
that results from the eigenevolution model has a slightly smaller a, especially for m < 0.5 Mq (thick solid 
histogram). This effect is larger for the default case {a^ — 2.3), because the larger number of massive stars 
implies more systems in which the secondary gains mass as a result of eigenevolution. The effect on a is too 
small, however, to make a significant difference in the alpha-plot (e.g. Fig. g below). Fig. ^ also displays the 
large difference between the system MF and the single-star MF at low masses. The IMF has a maximum 
near 0.1 Mq, whereas the system MF has one near OAMq, and underestimates the number of 'stars' by an 
order of magnitude near m = 0.01 Mq, and by a factor of three near m = 0.08 Mq. 



4.4. Nota bene 

The cluster models constructed here are extremes, in that they have a very high central density equal 
to that observed in the ONC. This assumption leads to a rapid depletion of the binary population, as 
shown below (Fig. M; see also de La Fuente Marcos 1997). Disruption of binaries occurs on a crossing-time 
scale (Kroupa 2000a) in any cluster, so that it takes much longer in real time for the binary population 
to decrease in a Pleiades-type cluster, for which Kahler (1999) shows that / « 0.7 is possible. Likewise, 
the pre-main-sequence cluster IC348, which has a density of about 500 stars/pc'^, has a binary proportion 
similar to that in the Galactic field (Duchene, Bouvier & Simon 1999). As shown by KTG91, such a binary 
proportion requires significant correction to the observed system LF to infer the IMF. The problem with 
unresolved binaries may be even worse for lower-density clusters still, such as studied by Testi, Palla & 
Natta (1999), because the binary population evolves on much longer time-scales, and is thus likely to be less 
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evolved than in the clusters studied here. The problem will never be smaller in such clusters, unless they 
consist of a stellar population that had an unusually small initial binary proportion (/ < 0.3), i.e. smaller 
than even in the evolved extreme models here. Such a population has never been observed in any Galactic 
cluster or association to this date (e.g. Ghez et al. 1997; Duchene 1999). 

Any real population is thus likely to have a larger binary proportion than in the models considered here 
after about three crossing times (~ 0.8 Myr). In addition, the present results will be an underestimate of 
the bias, because only binary systems are considered. Real populations contain something like 20 per cent 
or more triple and quadruple systems, which, when not resolved, increase the systematic error made in the 
observational estimate of a. What is inferred in this paper is thus the minimum correction to a. 

This is particularly true for to<; 1 A/q, because the observed mass-ratio distribution for massive stars 
(e.g. Preibisch et al. 1999) has secondaries that are typically more massive than IMq, whereas in the 
models here, massive primaries typically have very low-mass companions owing to the random sampling 
hypothesis. This is very important when considering the system MFs below. It will be evident that the 
models lead to essentially no bias for massive stars, but this is more likely to be a shortcoming of the present 
assumptions, rather than proving that the IMF for massive systems is not subject to a significant bias, as 
discussed in Section p.l.l| . Clearly, this is a fundamentally important topic requiring much more work to 
construct a more realistic initial mass-ratio distribution for massive stars. In addition, a systematically 
different IMF between the LMC and the MW for massive stars may become evident, if the binary properties 
differ systematically between the two galaxies, because then the correction for systematic bias would be 
different for the two samples. At present no such difference is known, and so the empirical LMC and MW 
data plotted in Fig. |l| can, at present, be only taken to mean that the IMF for massive stars may be the 
same in the two populations. 



5. RESULTS 

The results obtained from the stellar-dynamical calculations are used to study temporal and spatial 
apparent variations of the single-star and system MFs. 



5.1. Cluster evolution 

As an impression of the evolution of the star clusters. Fig. H displays the scaled number of systems and 
single stars with R < 3.2 pc. iVsys(t) = iVsing(i) + N\)in{t) increases for t^ 2.5 Myr because the disruption of 
binary systems liberates secondaries. That is, the observer would find that the number of 'stars' increases 
with time. After t w 2.5 Myr, Ngys decreases with a rate depending on N, because the clusters expand 
owing to binary-star heating, relaxation and mass- loss from evolving stars. 

The binary proportion (Fig. pi) decreases within a few initial crossing times. The decay occurs on 
exactly the same time-scale for the different clusters, demonstrating that it is not the velocity dispersion in 
the cluster alone which dictates the disruptions, but the density as well. Owing to the ejection from the 
cluster of preferably single stars and because of mass segregation, / is larger for systems with R < 3.2 pc 
and at times t^2.5 Myr, than for systems at larger distances from the clusters. The least massive clusters 
(N = 800) have expanded appreciably by this time so that the remaining binary population in the cluster 
is hard, and no further significant disruption of binaries occurs (/ ~ 0.55 and increasing for t^2.5 Myr). 
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The more massive clusters, however, remam more concentrated for a longer time (top panel of Fig. H), and 
consequently the binary-star hard/soft boundary remains at a higher binary binding energy for a longer 
time. At any time t^Sicross ~ 0.7 Myr, the binary proportion is higher in the clusters with smaller iV, 
which his particularly evident in Fig. |6| below. This is a nice example of the caveat raised in Section 4.4. 



Further details on these processes are available in Kroupa (2000b), and in the seminal paper by Heggie 
(1975). 

The evolution of the binary proportion for primaries with different masses is illustrated in Fig. |6[ 
The binary proportion of BDs falls rapidly, and stabilises near /bd = 0.20 in all models. M dwarfs retain 
a much higher binary proportion by t = 150 Myr, /m ~ 0.4 — 0.5, depending on N, and more massive 
primaries retain a slightly higher binary proportion still. The overall binary proportion of O primaries 
(m > 8Mq) shows a complex behaviour. Initially, most O primaries have low-mass companions. These are, 
however, exchanged for more massive companions near the cluster core. When the primaries explode, these 
companions are left or are ejected as single stars. In addition, violent dynamical encounters in the cluster 
core eject single massive stars. Overall, /o decays, but higher-order multiplicities that form in three-body 
encounters are not accounted for. 



5.2. The alpha plot for cluster populations 

Having briefly discussed the evolution of the clusters and of the binary population, the following question 
can now be posed: What MFs would an observer deduce if an ensemble of such clusters were observed 
at different times, under the extreme assumption that the mass of each star or system can be measured 
exactly? 

Figs. to H show the results for each N. The upper panels assume the observer sees all stars with 
R < 3.2 pc, whereas in the lower panel it is assumed that only the system masses can be measured exactly 
for systems with R < 3.2 pc. The MFs are constructed at times t = 0, 3 Myr and 70 Myr. For the 
single-star MFs, the results at t = are the same as for pure statistical noise (Fig. ^). 

At i = 0, the single-star IMF is well reproduced. The system MF, on the other hand, underestimates 
a significantly for nisys < 1 Mq, with oq « —0.8 (instead of H-0.3) for to ^0.08 A/0, ai « +0.7 (instead of 
+1.3) for O.O8<TO<O.5Af0, and a 2 w +1.5 (instead of +2.3) for 0.5<to<1Mq. 

At t = 3 Myr and 70 Myr, most of the BD systems have been disrupted (Fig. ||), with typically 
/bd ~ 0.2, and most star-BD systems have also ceased to exist, so that ao is only slightly underestimated 
for the system MF. Work is in progress to study if the resulting mass-ratio distribution becomes consistent 
with the observed 'BD-companion dessert' for nearby stars (M. Mayor 2000, priv.commun.). In mass ranges 
&2 and 63, the power-law index is still underestimated significantly, because the surviving binary proportion 
is typically / > 0.4 for to > 0.08 Mq. For 62 the lower panels in Figs. to || read ai ~ +0.8, and in 
63, a2 ~ +1.7. The bias in measuring ai,2 for the system MF rather than the single-star MF is thus not 
significantly reduced at later times. 

This bias will operate for even older clusters, because further binary disruption is essentially halted 
in the expanded clusters, and / begins to increase with time as energy equipartition retains the heavier 
binaries in the cluster at the expense of single stars (fig. 3 in Kroupa 1995d). However, with time the bias 
will decrease for 02 as the turnoff mass becomes smaller, i.e. as the number of primaries with to ^ 1 Mq 
decreases. As an extreme example, globular clusters retain a significant proportion of their low-mass 
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stars (Vesperini & Heggie 1997), but stars with m^Q.8 Mq have ceased to exist, so that no m'^0.8 Mq 
companions are 'hidden' with brighter primaries. 

For N = 800 (Fig. H) the scatter in range 65 is very large, and rather similar to what is seen in the 
observational data in the shaded area (0.8 — SMq, Fig. |l|). This is interesting, because in these models it is 
the stars in the mass range 3 — 9 Mq that are the most massive and abundant enough to eject each other 
from the core after meeting there through mass segregation, causing large fluctuations in the measured MF. 
The same holds true for the cluster data in the shaded range in Fig. |^. For example, p Oph contains not 
more than a few hundred systems, so that the most massive stars populate roughly the shaded range. The 
Pleiades is 100 Myr old, so that stars with m^ 10 Mq have evolved off the main sequence, and the stars 
just below this mass interact in the cluster core. 

In summary, comparison of the three figures shows that the scatter in a decreases as N increases, but 
that the scatter is larger than pure Poisson noise (cf. the t > data in the upper panel of Fig. M with the 
N = 10"* model in Fig. 0). The most important result though is that ai_2 is underestimated by Aa ~ 0.5 
for the system MF in the mass range 0.1 — IMq. And, an observer deduces fewer BDs in an unevolved 
population (t = 0; Fig. P such as in Taurus-Auriga, than in a population that is older than a few crossing 
times, such as the Trapezium and the Pleiades Clusters (see also Kroupa, Aarseth & Hurley 2001). The 
figures also show that for a single-age population the scatter is always smaller for to ^ 1 Mq. For to ^ 1 Mq 
the scatter for the clusters with A^ = 3000 and A^ = 10'' stars is comparable to the observed scatter. Even 
when A^ — 10^, models with a^ — 2.3 cannot be differentiated from models with a^ — 2.7 in mass ranges 55 
and 66. 



5.3. The alpha plot for cluster halo populations 

The MF "in" young rich clusters can often only be determined by avoiding the crowded central regions. 
This can cause systematic uncertainties because stellar encounters lead to preferentially lower-mass stars 
and preferentially single stars populating an extended halo, or being ejected from the cluster. 

The clusters with A^ = 3000 and A^ = 10** stars are used to investigate the MF for systems lying at a 
distance R > 3.2 pc from the cluster centre. The results are shown in Fig. 11^, assuming the observer can 
only determine the masses of systems. 



The scatter is larger than within the clusters {R < 3.2 pc. Section |5.2| ), and the bias for m < 0.5 Mq 
that leads to an underestimate of ai in binary-rich populations is reduced significantly. This results because 
the halo population is depleted in binary stars (Fig. ph. 

Two extreme examples are marked with double symbols. The corresponding MFs are plotted in 
Fig. [l^. One example is the system MF for a halo population at an age of 3 Myr. It's particularly flat MF 
for ?7i > 10 Mq {a = 0.97) comes about because the cluster core just expelled a few massive stars to the 
outer regions. The steep MF for a 70 Myr old population with a = 4.85 at iTOav = 0.9 (double square in 
the lower panel) arises because stellar evolution has removed stars with tti^ 10 Mq, and because the stars 
with a mass just below the turn-off mass are located preferably near the cluster core. The fitted power-law 
indices are listed in Table g. 
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logio(m) [Mo] 



Fig. 4. — Mass functions for single-stars (solid histograms) and systems (dash-dotted histograms) at i = 
in models B1E4 (thick histograms) and BlE4d (thin histograms). Note the smaller number of massive stars in 
model BlE4d, which has a steeper IMF for m > IMq with as = 2.7 (Table |). The solid dots are the IMF for 



N = 10' 



stars (Fig. g) scaled to iV = 10'', and the vertical dotted lines and arrowed regions are as in Fig. ^ 
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Table 3: The two examples highlighted in Fig. hol The corresponding MFs are plotted in Fig. hll. The table lists 
the number of log- mass bins used in the fit, nb, the average log-mass over which the fit is obtained, logjQTTiav, the 
fitted power-law index a and the probable uncertainty (Tf_a. 
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Fig. 5. — Examples of the evolution of individual clusters. Top panel: The number of systems (thick curves) and 
all individual stars and BDs (thin curves) within the innermost 3.2 pc. The short-dashed lines are for A'^ = 800, the 
solid lines are for A'" = 3000, and the dash-dotted lines are for A*' — 10'*. Bottom panel: The binary proportion for 
R < 3.2 pc (thick curves), and all R (thin curves) for the same cases as in top panel. In both panels, the horizontal 
dotted lines indicate the times (3 and 70 Myr) at which the mass functions are observed. 
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Fig. 6. — The proportion of binaries with primary masses > 8Mq, /o (thick long-dash-dotted curve), 1 — 8Mq, /im 
(thin short-dash-dotted curve), and 0.5 — IMq, /k (solid curve). M dwarf primaries (0.08 — 0.5 Mq) have a binary 
proportion, /m (thin dashed line), whereas brown dwarfs (0.01 — 0.08 Mq), /bd, are shown as the thin dotted line. 
Each curve is an ensemble mean. 
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Fig. 7. — The Q-plot for five B800 models at t = 0, 3 Myr (left panels) and 70 Myr (right panels). The single-star (top panels) and system (lower 
panels) MPs are constructed for stars with R < 3.2 pc. The open circles, open triangles, vertical and horizontal lines are as in Fig. 3. 
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Fig. 8.— As Fig. 7, but for five B3000 models. 
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Fig. 9. — As Fig. 7, but for models with A'" = 10* stars. The crosses, four-pointed stars and squares are for the two B1E4 models, whereas the same 
symbols but rotated by 45 degrees are for the two BlE4d models. 
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R>3.2pc, System MF 
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Fig. 10. — Spatial variation of the MF: The MFs in models B3000 (5 renditions), B1E4 (2 renditions) and BlE4d 
(2 renditions) for systems with R > 3.2 pc. Two particularly exotic examples are highlighted using double symbols, 
and the corresponding MFs are plotted in Fig. nil, and the MF-fits are listed in Table pi Otherwise as Fig. bl Note 
the changed a-scale. 
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Fig. 11. — Spatial variation of the MF: The MPs for systems with R > 3.2 pc for model B1E4 showing two cases: 
solid histogram, f = 3 Myr (four-pointed double star in Fig. hd) and dash-dotted histogram, t — 70 Myr (double 
square in Fig. |lO[ ). The open and filled circles represent the A'^ — 10*' star IMF from Fig. H after appropriate scaling. 
The arrowed mass-ranges are as in Fig. m. 
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5.4. A synthetic alpha plot 

The results from all cluster models at different times and for the inner and outer cluster regions can be 
combined to form a synthetic ensemble of populations. The result is shown in Fig. 02 for the case that the 



observer is able to measure the mass of each star exactly. Fig. 13 shows the results assuming the observer 
can measure the system masses exactly. 

The model a values obtained by fitting power-laws to the model system MFs are consistent with 
a ~ 2.3 for to^IMq, thus re-deriving the input IMF despite unresolved binary systems. This result will 
be re- visited in future work for the reasons stressed in Sections |2.1.1| and 4.4. 



For m < lA^Q, the average system-a are too small, except in the BD regime, where approximately 
the input value (ao = +0.3) is arrived at because of the small surviving binary proportion. Fig. |l3| 
thus demonstrates that the observational data (open circles and triangles) underestimate the single-star 
power-law index in mass-ranges 62 and 63 (eq. 0) by about Aa w 0.5, because binary systems are not 



resolved. This is a reliable result, because of the reasoning in Section 4.4, that is, because the cluster library 
used here has an extreme initial density. Any Galactic embedded cluster with a lower density may lead to 
a larger bias, because in lower-density clusters the binary population is eroded at a slower rate, allowing a 
higher binary proportion to survive for longer times. The binary proportion is certainly not lower in such 
clusters, which is also confirmed by detailed analysis of observations (e.g. Kahler 1999 for the Pleiades; 
Kroupa & Tout 1992 for the Praesepe). 

Again, it is stressed that the above corrections to a are minimum values, especially for BDs. The binary 
proportion of these may be larger in clusters with lower density, because it takes longer for /bd to decrease 
in lower-density clusters. The maximum corrections to be applied to the observed, i.e. system MFs, are 
derived from the models at t = (e.g. Fig. |): Aa = +1.3 for BDs, and Aa = +0.8 for 0.08 < m < 1 Mq. 
Such large corrections are, however, unlikely, because / < 1 usually (except in Taurus-Auriga, cf. Luhman 
2000). 

The observational data in Fig. |l| therefore imply a single-star IMF that is steeper than eq. g for 
O.O8^m^lM0 by Aa ~ 0.5 at least. Thus, for these data the corrected IMF has ai « 1.8 for 
0.08 — 0.5 Mq, and a2 ~ 2.7 for 0.5 — 1 Mq, probably with unchanged ao and a^. The implications of this 



are discussed in Sections 6.3 and 6.3, 



Figs |l2| and |l^ show that the model scatter in a is similar to that seen in the observational sample. 
Despite starting in each case with the same IMF, an observer deduces power-law indices that have a scatter 
of about (Tq = 0.5 for m'^1 Mq and CTq = 1 for m ^ 1 Mq, even if each stellar mass is measured exactly. 
The finding is thus that the IMF can never be determined more accurately than this scatter, and that the 
scatter seen in the alpha-plot (Fig. |l|) can be explained with Poisson noise and stellar dynamical effects. 



6. DISCUSSION 

A cautionary remark concerning the alpha-plot is made, namely that in reality the left and right parts of it 
are disjoint. Also, some tentative evidence for a systematically varying IMF is presented, especially in view 
of the proposed revised IMF. 
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All models, Stellar MF, t = 3 and 70Myr 
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Fig. 12.— All models B800 (5 renditions), B3000 (5 renditions) and B1E4 (2 renditions) for i = 3 Myr and 70 Myr 
for individual stars with R < 3.2 pc and R > 3.2 pc. Only power-law fits that are based on more than nb = 3 log-mass 
bins are plotted. The horizontal and vertical lines, the faint open circles and open triangles have the same meaning 
as in Fig. H. 



All models. System MF, t = 3 and 70Myr 



a 2 




logn,m [Mq] 
Fig. 13. — As Fig. h3, but assuming the observer cannot resolve systems. 
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6.1. The dichotomy problem 

When considering the alpha-plot (Fig. |l|), it must be remembered that that left (jTi^lAf©) and right 
(to ^ 1 Mq) parts of it are actually disjoint. 

That is, any nearby cluster that is older than a few Myr so as to allow the application of reasonably 
well understood pre-main sequence or main sequence stellar models, contains no O stars or is already too 
old for them to still exist. This is very true for the Galactic-field IMF - there is only an indirect handle on 
TO ^ 1 Mq stars through stellar remnants, but this requires an excellent understanding of stellar evolution, 
the sfh and Galactic-disk structure (e.g. Scalo 1986). Conversely, any population of stars for which the MF 
is constrained through observations for to > 1 M© is usually so far away that the left part of the alpha-plot 
is not accessible to the observer, and/or so young that measuring the derivative (a), i.e. the shape, of the 
IMF for m < 1 Mq becomes a lottery game because of the uncertain pre-main sequence tracks (Section ||). 

That low-mass stars do form in large numbers in any population that also forms O stars is established. 
Examples are the ONC (Hillenbrand 1997), R136 in the 30 Dor region in the LMC (Siriani et al. 2000), 
and NGC 3603, the most massive visible HII region in the MW (Brandl et al. 1999). However, the ONC 
is so young that mass estimates become unreliable, compromising conclusions about the detailed shape of 
its IMF, and in the other cases the census of low-mass stars is not complete. Thus, the shape of the IMF 
spanning logioTO, — —2 to 2 is not known for any population, and it remains an act of faith to assume that 
the IMF can be approximated by the form of eq. ||. 

Globular clusters consist entirely of low-mass stars today, but the existence of neutron stars 
demonstrates that massive stars formed in them as well. Paresce & De Marchi (2000) suggest that the 
MF for a sample of a dozen globular clusters can be fit by a log-normal MF with approximately one 
characteristic stellar mass and standard deviation. A further analysis will show how the differences compare 
with the spread in a seen in Fig. H^. More interesting in the present context is that Piotto & Zoccali 
(1999), who use the same stellar models by Baraffe et al. (1997) as Paresce & De Marchi, demonstrate that 
power-law MFs fit rather well for a wide range of globular clusters, with a « 0.5 — 1.2 for m ^ 0.5 — 0.7 Mq, 
but the IMF is not measurable for stars with to,^0.7Mq. 



6.2. A revised IMF 



In Section 5.4 the suggestion is made that the systematic bias towards low ai and a2 due to unresolved 
binaries implies that the single-star IMF may be steeper than inferred from observations that do not resolve 
binary systems. Correcting the ensemble of observed a in Fig. |l| for this bias leads to the following revised 
IMF, 

ao = +0.3 ± 0.7 , 0.01 < to/Mq < 0.08, 

^, iiO-LAc: n r\Q ^ .^^ /]\/r ^ n r:n 

(6) 



ao 


= +0.3 ±0.7 


, 0.01 < tu/Mq < 0.08 


ai 


= +1.8 ±0.5 


, 0.08 < tti/Mq < 0.50 


0-2 


= +2.7 ±0.3 


, 0.50 < to/Mq < 1.00 


as 


= +2.3 ±0.7 


, 1.00 < tu/Mq, 



where the uncertainties from eq. are carried over. 

The revised IMF has, for stars with 0.01 < m < 50 Mq, an average stellar mass < to >= 0.20 M© and 
leads to the following population: 50 % BDs (0.01 - 0.08 M©) contributing 10 % to the stellar mass, 44 % 
M dwarfs (0.08 - 0.5 Mq) contributing 39 % mass, 4.3 % "K" dwarfs (0.5 - 1.0 Mq) contributing 14 % 
mass, 2.3 % "intermediate mass (IM) stars" (1.0 — 8.0 Mq) contributing 24 % mass, and 0.15 % "O" stars 
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(> 8 Mq) contributing 12 % mass. O and IM stars thus contribute together 36 per cent of the total mass. 
If q;4 = 1.15 (m > 8 Mq) then 50 per cent of the mass is in stars with 8 < ?ti < 120 Mq. 

This revised IMF can be viewed as the present-day star-formation IMF, and is in good agreement with 
the pre-stehar clump MF measured by Motte at al. (1998) and Johnstone et al. (2000) for p Oph: ai « 1.5 
and a2 ~ 2.5; especially so since each clump is likely to form a multiple star. 



6.3. Possible evidence for a variable IMF 

A short account is made of the most promising evidence for a systematically varying IMF. The discussion 
in Sections 6.3.1 to 6.3.2 is visualised in Fig. |lj, in which the various IMFs are compared. 



6.3.1. Globular clusters vs Galactic field 

The suggestion in Section ^^ that the alpha-plot (Fig. |^) may imply a present-day star-formation (pdsf) 
IMF (eq. g) that is steeper than the Galactic-field IMF (eq. g) is interesting when compared to the MFs 
estimated for globular clusters (Section ^J). These are very ancient and metal poor systems, so that 
a systematically different IMF (Larson 1998) ought to be manifest in the data. The difference should 
be in the sense that globular clusters ought to contain a characteristic stellar mass that is larger than 
that in more metal-rich populations. The systematically flatter MF in globular clusters compared to the 
Galactic-field IMF (eq. ||), and especially to the pdsf IMF (eq. ^, may thus be due to a real difference in 
the star- formation conditions. 

However, unfortunately the evidence is not conclusive because globular clusters have lost preferentially 
low-mass stars, leading to a systematic fiattening of the MF with time, unless the clusters are at large 
Galactocentric distances (Vesperini & Heggie 1997). The binary proportion in globular clusters is typically 
smaller (/^0.3) than in the Galactic field (/ « 0.6) but probably not negligible (Hut et al. 1992; Meylan & 
Heggie 1997), and correction for their effects may also steepen the measured MF. Approximate corrections 
that increase the measured a are Aa < 1 for dynamical evolution (fig. 6 in Vesperini & Heggie 1997) and 
Aa ~ 0.2 for unresolved binary systems, but a case-by-case study is required for detailed estimates. In their 
sample, Piotto & Zoccali (1999) find evidence for flatter MFs at smaller Galactocentric distances suggesting 
loss of low-mass stars as being an important bias. But, there is also evidence for a correlation such that 
more metal-rich clusters have larger a. 

The Galactic-field IMF (eq. ||) is valid for stars that are, on average, about 5 Gyr old, and which were 
formed at a different epoch of Galactic evolution than the stars in the clusters featuring in Fig. |l|. This, 
then, suggests a possible systematic shift of star formation towards producing relatively more low-mass 
stars as star-formation moves towards conditions that may favour lower fragmentation masses through 
higher metallicities and lower cloud temperatures. That the pre-stellar core MF in p Oph is somewhat 
steeper than the Galactic-field IMF (eq. ||), while being consistent with a fragmentation origin (Motte et al. 
1998), supports this notion. 
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A systematically varying IMF ? 



!2fl 
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logjom [Mq] 



Fig. 14. — Evidence for a systematically changing IMF. The present-day star-formation (pdsf ) IMF (eq. m) is shown 
as the thick dashed line. The Galactic-field IMF (eq. El) is the thick solid line. It is truncated at m = 1 Mq to 
express our ignorance about the IMF for m > 1 Mq for this population that has an average age of about 5 Gyr 
(the "dichotomy problem", Section 3.1). In both cases the shaded areas represent the approximate 95-99 per cent 



confidence region. For comparison, the Miller & Scalo (1979) log-normal IMF for a constant star- formation rate and 
a Galactic disk age of 12 Gyr is plotted as the thin long-dash-dotted curve (it's derivative is shown in Fig. hi). Seven 
globular clusters give a\ — 0.89, with upper and lower values of 1.22 and 0.53, and Q2 = 2.3 for 0.6 < m < 0.8 Mq 
(Piotto & Zoccali 1999) as indicated by the short-deished lines and the heavily shaded area. Three possible IMFs 
for Galactic-halo WD-progenitors are suggested by the thick long-dash-dotted and short-dash-dotted lines (Chabrier, 
Segretain & Mera 1996, CSM96), and the thick dotted line (Adams & Laughlin 1996, AL96). The MFs have been 
scaled such that they agree near 0.5 M0, except for the ancient IMFs, which are scaled to fit the Galactic-field IMF 



near 2 M, 



©• 
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6.3.2. Galactic-halo white dwarfs 

Another possible empirical hint for a variable IMF may be provided if part of the dark halo of the Galaxy 
were in the form of ancient white dwarfs. This is becoming a distinct possibility, given that a handful of 
candidate ancient halo white dwarfs have been discovered (Elson, Santiago & Gilmore 1996; Ibata et al. 
1999, 2000; Mendez & Minniti 2000). 

From eq. g one obtains per WD progenitor (1 — 8 A/q, e.g. Weidemann 1990) about 8 dwarfs with 
TO — 0.1 — 0.7 Mq. No such halo dwarfs that might belong to the same population as the putative WDs 
have been found, requiring a radically different IMF for their progenitor stars than is seen today in the 
Galactic disk. Also, for consistency with chemical enrichment data, such an IMF cannot have many stars 
with TO>5Afo (Chabrier, Segretain & Mera 1996; Adams & Laughlin 1996; Larson 1998; Chabrier 1999). 



6.3.3. Radial variation in a very young cluster 

Hillenbrand (1997) demonstrates that the ONC has pronounced mass segregation, and this may be 
interpreted as an IMF which has a radial variation, if dynamical mass segregation is not fast enough to 
produce such mass segregation within the age of the cluster. The age of the ONC is estimated to be less 
than 1 Myr for most ONC stars (Hillenbrand 1997; Palla & Stabler 1999), and Bonnell & Davies (1998) 
suggest, by using a softened A^-body code, that mass segregation takes too long to produce the observed 
effect. However, stellar-dynamical computations with a direct A^-body code that correctly treats the 
many close encounters must be applied to this problem (Kroupa, in preparation). If the mass-segregation 
time-scale is too long to produce the observed effect, then we would have a well-documented case of a 
variable IMF most likely through interactions of pre-stellar cores, as suggested by Bonnell et al. (1998) and 
Klessen (2001). 



7. CONCLUSIONS 

The following three main points are covered in this paper: 

I. The Galactic-field IMF. The form of the average IMF consistent with constraints from local star-count 
data and Scalo's (1998) compilation of MF power-law indices for young clusters and OB associations is 
inferred. The IMF is given by eq. g. This form may be taken as the universally valid IMF. 

II. The alpha-plot: scatter and systematics. Assuming the universal IMF (eq. H), how large are the apparent 
variations produced by Poisson noise, the dynamical evolution of star-clusters and unresolved binary 
systems? 

This is studied by making use of the alpha-plot, in which IMF power- law indices inferred for A^-body 
model populations are plotted as a function of stellar mass. The extreme assumption is made that the 
observer can measure each stellar or binary-system mass exactly. The resultant apparent variation of the 
IMF thus defines the fundamental limit for detecting true variations. Any true variation of the IMF that 
is smaller than this fundamental limit cannot be detected. This is the reason why no robust evidence for a 
variable IMF has surfaced to date. The available population samples are too small (e.g. one ONC vs one 
p Oph). 

The model clusters have an initial binary proportion of unity and contain A^ = 800, 3000 and 10"* stars 
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with a central density as in the ONC. Clusters with a smaller initial density evolve on a longer time-scale. 
The binary-star problem is thus potentially worse in less-dens clusters, because binary systems survive for 
longer. 

The observed spread of power-law indices is arrived at approximately. For the ensemble of model 
clusters studied here it is (7^ « 0.7 for BDs, aa ~ 0.5 for stars in the mass range 0.1 — 1 Mq, and ctq « 1 for 
stars with m > 1 Af© (Fig. |l|). 

For stars with to ^ 1 Mq, the system MF has, on average, the same power-law index as the underlying 
single-star IMF. That is, the present models do not lead to any systematic bias in this mass range (but see 
caveat in Section ^.l.l| ). Similarly, for BDs the input a is arrived at in the mean, but only if the population 
is at least a few crossing times old, because by then most BD binaries and star-BD binaries have been 
disrupted. For a dynamically younger population, a (and the number of BDs) will be underestimated 
depending on the binary proportion. 

To correct for unresolved binaries, the measured power-law index has to be increased by 0^ Aao £ 1.3 
for BDs and 0.5^ Aai^2 ^0.8 for 0.08 < to < IMq, the upper and lower Hmits applying for clusters 
that are unevolved (i = 0) and a few crossing times old, respectively, assuming / = 1 when t = 0. For a 
population in a cluster that is a few crossing times old, the corrections reduce to Aao « and Aai.2 ~ 0.5. 
These corrections have to be applied to any young population to infer the single-star IMF. 

Finally, as a cautionary remark, the left and right parts of the alpha-plot are observationally disjoint. 
It is an act of faith to assume that a{m) has the smooth dependence given by eq. y. 

III. IMF variations. Applying the above corrections to the ensemble of observed young clusters, a revised 
(or present-day star- formation) IMF is arrived at (eq. |g) . It is steeper for to ^ 1 M© than the Galactic-field 
IMF (eq. g), which is a mixture of star- formation events with an average age of about 5 Gyr. The pre-stellar 
clump mass-spectrum in the present-day star-forming cloud p Oph (Motte et al. 1998; Johnstone et al. 
2000) also indicates a steeper single-star MF than the Galactic-field MF. Intriguingly, the ancient MFs in 
globular clusters have a^O, but closer to than the Galactic-field IMF. The recent detection of candidate 
white dwarfs in the Galactic halo suggests that the IMF of the progenitor population must have been 
radically different by producing few if any low-mass and massive stars (a << for m'S^O.5 Mq and a >> 
for to>2Mq). 

Furthermore, the well-developed mass segregation in the very young ( ^ 2 Myr) ONC may exemplify a 
locally radially- varying IMF, if dynamical mass segregation is too slow. If A^-body calculations confirm this 
to be the case (work is in progress), then the ONC will be definite proof that the local conditions determine 
the average stellar mass, rather than it merely being the result of statistical fluctuations. 

The tentative suggestion is thus that some systematic variation may have been detected, with 
star-formation possibly producing relatively more low-mass stars at later Galactic epochs. Such a variation 
would be expected in the mass range ('^IMq) in which turbulent fragmentation, which depends on the 
cooling rate and thus metal abundance, dominates. Future observations of LMC populations might verify 
if the IMF has systematically smaller a for to ^ 1 Mq than the Galactic-field or present-day star-formation 
IMF. Unfortunately though, even if there is a trend with metallicity, it will be very arduous to uncover a 
systematic difference in a between the MW and LMC at low masses because the metallicity difference is 
not very large while the a-scatter is. A lack of systematic differences in a for to, ^ 10 Mq between MW 
and LMC populations may be a result of one physical mechanism, such as coalescence, dominating in the 
assembly of massive stars (Larson 1999). 
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